Spin-orbit coupling and intrinsic spin mixing in quantum dots 
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Spin-orbit coupling effects are studied in quantum dots in InSb, a narrow-gap material. Compe- 
tition between diflterent Rashba and Dresselhaus terms is shown to produce wholesale changes in 
the spectrum. The large (and negative) g-factor and the Rashba field produce states where spin is 
no longer a good quantum number and intrinsic flips occur at moderate magnetic flelds. For dots 
with two electrons, a singlet-triplet mixing occurs in the ground state, with observable signatures 
in intraband FIR absorption, and possible importance in quantum computation. 
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The creation and manipulation of spin populations in 
semiconductors has received a great deal of attention in 
recent years. Conceptual developments that have mo- 
tivated these efforts include prominently the Datta-Das 
proposal for a spin field-effect transistor,^ based on the 
Rashba spin-orbit coupling of electrons in a 2DEG)^ and 
the possibility of building quantum computation devices 
using quantum dots (QDs) J?^ It is then important for full 
control of spin-flip mechanisms in nanostructures that all 
spin-orbit (SO) efi^ects be understood. 

There are two main SO contributions in zincblende 
materials like in addition to the structure inver- 

sion asymmetry (SIA) caused by the 2D confinement 
(the Rashba effect), there is also a SO term caused by 
the bulk inversion asymmetry (BIA) in those structures 
(the Dresselhaus term)ii^ Notice that additional lateral 
confinement defining a QD introduces another SIA term 
with important consequences, as we will see in detail. 
Although the relative importance of these two effects de- 
pends on the materials and structure design (via intcrfa- 
cial fields) , only recently have authors begun to consider 
the behavior of spins under the influence of all effects. 
For example, a modification of the Datta-Das design was 
recently suggested to allow for a diffusive version of the 
spin FET)4 and that proposal relies on the additional 
influence of the Dresselhaus SO coupling in the system. 

Work in wide-gap materials (mainly GaAs))& uses 
a unitary transformation on the Hamiltonian of the 
system,- after which one gets an effective diagonal SO 
term which incorporates the Rashba effect in a pertur- 
bative fashion. That approach is valid since the SO cou- 
pling is small in GaAs. However, the approximation is 
not valid for all of the A-^B^ structures, as it is the case 
for InSb, for example, where both SIA and BIA effects 
are anticipated to be much larger.^ In this material, one 
needs to deal with the full Hamiltonian. 

There are just a few works discussing SO effects in 
narrow-gap nanostructures. Among them, uses k • p 
theory in InSb QDs in order to include SIA SO terms 
from both the Rashba field and the lateral confinement 
which defines the QD. This last SIA term is considered 



in and since it is diagonal in the Fock-Darwin (FD) 
basis no level mixing is found nor expected. In contrast, 
level mixing events are clearly identified in 0. Exper- 
iments in InSb QDs have explored the FIR response in 
lithographically defined dots,^*' and PL features of self- 
assembled dotsiii 

The goal of this work is to show how important differ- 
ent types of SO couplings are in the spectra of parabolic 
QDs built in narrow-gap materials such as InSb. We 
consider the Rashba-SIA diagonal and SIA non-diagonal, 
as well as the Dresselhaus-BIA terms in the Hamilto- 
nian, and proceed with its full diagonalization, in order 
to study features of the spectrum as function of magnetic 
field, dot size, (/-factor, and electron-electron interaction. 
We draw attention to the appearance of strong level an- 
ticrossings (mixing) for moderate magnetic fields in typ- 
ical QDs, and how this phenomenon (and 'critical' field 
where it occurs) is modified by the BIA terms not con- 
sidered before.^ As the level mixing involves states with 
different spin, this induces strong intrinsic spin flips in 
the system, regardless of the strength of the SO coupling, 
providing an important channel for spin decoherence in 
these systems. Moreover, measurement of FIR absorp- 
tion would yield direct access to the coupling constants; 
i.e., the dispersion of FIR absorption peaks and appear- 
ance of additional/split-off features are a direct conse- 
quence of the level mixing introduced by SO. 

Model. Assuming a heterojunction or quantum well 
confinement V{z) such that only the lowest z-subband 
is occupied, the Hamiltonian in the absence of SO inter- 
actions for a QD further defined by a lateral parabolic 
confinement is given by Hq = ^^k^ + V(p) -\- ^g^B^ ■ cr, 
where k — — iV + eA/{hc), and the in-plane vector 
potential A = i?p(— sin0, cos6', 0)/2 describes a per- 
pendicular magnetic field B = _Bz; m is the effective 
mass in the conduction band,^^ g is the bulk g-factor, 
/is is Bohr's magneton, V{p) = ^mojQp^ is the lat- 
eral confinement with frequency ujq, and in the Zee- 
man term ax,Y,z are the Pauli matrices. The analyti- 
cal solution of Hq yields the FD spectrum with energies 
Enia = (2n-|-|;| + l)fif7+Zftwc/2+.g/iB-Bcr/2, with effective 
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frequency = + ^^^74, and cyclotron frequency 

ujc = eB/{mc); states are given in terms of associ- 
ated Laguerre polynomials^ The confinement, magnetic 
and effective lengths are, respectively, Iq = ^ h/{muj^)), 

Ib — \J'hl [mujc) and A ~ ^Jh/ (mQ,). 

The SIA SO term with coupling parameter a is 
HsiA = acr • (Vy X k), where the total confinement 
potential is V{y) = V{p) + V{z). One can then write 
HsiA = Hfi + Hgjj^, where the diagonal contribution 
coming from the lateral confinement in cylindrical coor- 
dinates is H^'jj^ = a^az [Lz + j^^), with the adi- 

mensional radial coordinate x = p/A, and Lz = —id/d9. 

The Rashba term coming from the perpendicular con- 
finement field dV/dz is 

dV , 

Hr = -oij^ [a+L^A^ + (T-L+A+\ , (1) 

where L± = exp(±i^?), a± = (tTx±*CTy)/2, and operators 
A± = Td/dx + Lz/x + a;AV(2Z|). 

In zincblende structures one should also consider 
the BIA SO bulk Hamiltoniam^ After averaging 
in the 2;-direction, due to quantization, one gets 
Hbia = 7 [f^xkxky - aykykl] + 7 (k"^) [ayky - a^kx] + 
70'z {kz) (fe^ — ky), where 7 is the coupling parameter, 
the resulting first (second) term is cubic (linear) in the 
in-plane momentum, and the last term is zero because 
(kz) = 0; also, (fc^) ~ (tt/zq)^, where zq is the z-direction 
confinement length. One may write the BIA SO term as 
Hbia = H^ + H£, where the linear Dresselhaus contri- 
bution is given by 

H^^i^^[cT+L+A+-a^L^A^], (2) 

while the cubic contribution H£ can be expressed in 
terms of cr^Lj_ and a±L±, and different powers in x, 
d/dx, and Lz-^'^ Notice that under a finite magnetic field, 
the matrix elements with a±L± in H£ are not hermitian, 
and one needs to symmetrize them^i if the field is zero, 
this problem does not occur Ji^ 

For the electron-electron interaction H^e , an expansion 
in Bessel functions for |ri — r2|~^ is employed. The 
basis states are properly antisymmetrized, describing the 
unperturbed spin eigenstates. 

The general form of the various SO terms in the Hamil- 
tonian exhibit already interesting characteristics. For ex- 
ample, the magnetic field plays a role via its linear depen- 
dence in Hgjj^, Hji, and H^, or its B to B^ dependence 
in Hj^M. Most interestingly, this form of the Hamiltonian 
yields selection rules explicitly, dictating which levels will 
be influenced by the SO effects. For example, at zero 
field the diagonal SIA term splits the levels according 
to the total angular momentum j. The Rashba term in- 
duces a set of anticrossings in the FD spectrum whenever 
Al = ±1 — — A(T at finite field (due to the cr±iip terms; 
mostly negative Vs are affected since their magnetic dis- 
persions allow for crossings); the lowest anticrossing is 



between {n,l,a} = {0,0,—} and {0,-1,+}. The cubic 
BIA terms (with CipLj.) induce a set of anticrossings, 
which obey Al = =f3 and Act — ±1; the first one at 
low S-field involves the states {0,1,—} and {0,-2,-1-}. 
Terms with a±L± in H^ and H£ do not induce anti- 
crossings, but rather split and shift the spectrum due to 
matrix elements with = ±1 = Act. Notice that the 
matrix elements between states with different n's are in 
general non-zero, so that the full diagonalization involves 
mixings with various n- values. 

Results. The sequence of FD states of Hq starts at 
zero S-field with {n,l,a} — {0,0, ±}, followed by the 
degenerate set of {0, — 1,±} and {0,1, ±}i^ Spin and 
orbital degeneracies are broken by B and the states with 
negative / and positive ct acquire lower energies because 
of the negative g-f actor. The lowest energy level crossing 
is between states {0,0,—} and {0,-1,-1-}, and the field 
where it occurs in the FD spectrum is 

^° (3) 
MB ^m\g\{m\g\ + 2) 

where rh = m/mo. The moderate value of B^ is a direct 
consequence of the large \g\ in InSb.^^ For GaAs {\g\ — 
0.44, rh = 0.067), for example, this level crossing appears 
only at ~ 9.4T for a much smaller confinement, 

hujQ = 2meV, and in the region where Landau levels are 
well defined. Weaker confinement (smaller luq) shifts this 
crossing to lower fields. Notice that for 5 < 0, Hr mixes 
these states (and H^ shifts the crossing to higher fields). 
For 5 > it is H^ that would produce relatively stronger 
level anticrossings (and Hr would only shift the spectrum 
weakly), and it would then be absent in non-zincblende 
materials like silicon. 

The energy spectrum for InSb QDs with typical 
characteristicSfi^ and for the full Hamiltonian is pre- 
sented in Fig. lA vs. B field. The spectrum is obtained 
by direct diagonalization using a FD basis with n < 4 
(or ten energy 'shells'), i.e. 110 basis states. We have 
studied the progressive changes to the FD levels when 
including different SO terms in H . 

The diagonal Hgj^^ term shifts energies but does not 
change appreciably the position of the first crossing 
shown by the dotted lines in Fig. lA at ~ 2.6T. The en- 
ergy shifts induce two new crossings at low fields (inset 
d), since the SO orders states according to their total 
angular momentum j = I + s; the highest (lowest) state 
at zero field has j — 3/2 (1/2) in the second shell. At 
about 0.2T one recovers the 'normal' sequence of states: 
{0, -1, +}, {0,-1, -}, {0, 1, +}, {0, 1, -}. This competi- 
tion between SO and magnetic field is similar to the Zee- 
man and Paschen-Back regimes in atoms. We should 
note that this level ordering is observed in |12| . 

The non-diagonal Rashba contribution Hr introduces 
strong state mixing for any value of the a parameter 
whenever FD levels with Al = —Act — ±1 cross. This 
mixing converts the crossings at B^ to clear anticross- 
ings. Higher levels which satisfy these selection rules also 
anticross at nearly the same field. The field- width and 
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FIG. 1: A. Full Hamiltonian H spectrum vs. B field for InSb 
QD as in H. Highlights in dashed boxes: a shows zero-field 
splitting in second shell (affected mostly by BIA SO terms), 
and crossing at about 0.3T; compare with inset d with only 
SIA terms and two crossings at 0.02 and O.OGT, and much 
smaller zero-field splitting. Second crossing for this shell in 
full H is at 3.4r (e box), b and c indicate anticrossings (AC) 
induced by Rashba term with AZ = —Act = ±1; first AC 
(arrow) in b involves states {0, 0, — } and {0, —1, -|-} ({0, 1, — } 
and {1, 0, -I-} in c). Dotted lines indicate FD levels crossing at 
2.6r. B. Lateral size dependence. Dotted lines: SO zero-field 
splitting in a box on left panel. Solid lines: Be field of first 
AC in b box on left panel; inset shows splitting at that AC. 
Arrows at 190A show QD size for spectrum in A. H\ curves 
(squares) use O : H2 (triangles) and H-j, (circles) use same 
parameters but four times stronger Rashba field {H2) or twice 
as large zq (Hs). Both cases increase relative strength of SIA 
terms. Solid line with no symbol shows B^ in 

energy-amplitude (or level splitting) of the mixing is dic- 
tated by the value of a, while the value of Be where the 
anticrossing occurs is nearly unaffected by a. 

The cubic Dresselhaus contribution H£ induces anti- 
crossings (via CTzpLj.) and zero- field splittings {ct±L±) in 
the FD spectrum. The splittings are much smaller than 
those induced by the Rashba term and practically unno- 
ticeable in the spectrum, reflecting the smallness of the 
Ej^ for these parametersi^^ For the linear BIA contri- 
bution, however, the <7±L± terms in have a much 
bigger impact on the zero-field splittings, which can in 
principle be 'tuned' by changing the effective z-size, zq- 
alone induces such a strong mixing at low fields that 
one cannot identify the two Zeeman and Paschen-Back 
regimes. 

Notice in the full spectrum of H (Fig. lA) that the 
first group of anticrossings (for n — levels) induced by 
Hn is shifted to higher field due mostly to so that 
Bq Be — 3.3T (box b and lower arrow). The set of 
anticrossings at ~ 5.5T is also due to Hn and arises from 
the n = 1 level manifold (box c and upper arrow) . At low 
field, only a single crossing in the second shell at ~ 0.3T is 



FIG. 2: A: Spin z-component vs. B for lowest 31 states; b and 
c labels refer to boxes in Fig. lA. Higher energy anticrossings 
in each set are shifted to lower fields. If only SIA terms are 
considered (panel B), all spin mixing occurs at field B ~ 2.6T. 
C: Sz for lowest 7 states; full spin mixing at anticrossing. D: 
dot with Rashba field 4 times stronger {Hi and H2 defined as 
in Fig. 1). Increasing SIA SO produces stronger mixing. 

present and dominated by (box a; compare with inset 
d and notice second crossing in box e). The sequence of 
the first excited levels at zero field is j — 3/2 (1/2) for 
higher (lower) energy, while at higher energies both SIA 
and BIA terms cooperate to produce anticrossings (not 
visible at the resolution in Fig. lA). 

Figure 2 illustrates the importance of the level anti- 
crossings on the spin, as the expectation value of for 
each state is plotted vs. B. Figs. 2A and B include all 
states with E < 80meV (for full SO and only SIA terms, 
respectively), while figs. 2C and D focus only on the low- 
est seven levels. Although a large majority of states have 
Sz close to ±1/2, as one expects for pure states, there are 
significant deviations. The various SO terms mix levels 
close to accidental degeneracy points in the FD spectrum 
and produce the large deviations seen in the figure. 2C 
shows how Hn produces an intrinsic (i.e., no phonon- 
assisted) total collapse of the spin number for the low 
energy states in the QD. Although the ground state is 
nearly pure {Sz — 1/2, and more so at higher B), the first 
few excited states totally mix at Be — 3.3T. 2D shows 
how a stronger Rashba field {dV/dz = — 2 x lO^^eV/A) 
greatly widens the mixing region and lowers Be — 2.8T. 

One can further appreciate the intricate balance of SO 
terms under a magnetic field. We analyze how various 
quantities are affected by changes in the the lateral and 
vertical sizes, Iq and zq, or the Rashba field dV/dz, as 
shown on Fig. IB. The zero-field splitting (dotted lines) 
is dominated by the linear BIA contribution for any value 
of Iq here. Increasing zq strongly reduces the splittings 
because the Dresselhaus contribution weakens; the reduc- 
tion is even more drastic if one increases dV/dz, which 
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makes the Hn contribution bigger and can then cancel or 
suppress better the sphtting produced by Hgj^^. Some 
authors have considered the possibihty of tuning such 
SO terms to produce total cancellation of the zero-field 
splitting, although considering only Hj^ and H^.* How- 
ever, one also has to take into account Hgjj^ and H£ 
contributions, which may be important (the zero-field 
cancellation occurs at values of zq or Rashba field about 
ten percent smaller than with only the former terms). 
One should notice, in any event, that this change in pa- 
rameters only eliminates the zero-field splitting but not 
the anticrossing at finite field, and measurement of both 
quantities on the same sample could yield information on 
the relative strength of the a and 7 parameters. 

The anticrossing field Be (solid lines/symbols) de- 
creases with QD size, roughly according to lj3Jl, Bq ~ 
luq ~ I/^/Iq. a finite a slightly increases Be, but 
the BIA contribution considerably upshifts it, as men- 
tioned above. Increasing zq or dV/dz decreases Be- At 
lo = 320A {hiJo = 5meV), Be = 1.6T, while it shifts 
to 1.15T if dV/dz is four times larger or to 0.85T if zq 
is doubled, both cases decreasing the BIA contribution. 
These values are comparable to those in Q without in- 
cluding BIA terms (adjusting for differences in system 
parameter values). Anticrossings at such low fields may 
be interesting for applications due to easier access. 

The energy splitting at Be (inset in Fig. IB), has main 
contribution from the Rashba term for any dot size con- 
sidered, but the BIA reduces the splitting substantially. 
If Zq is changed from 40 to 8OA the splitting is enhanced 
slightly, but larger zq produces no significant changes. 
However, the splitting is drastically enhanced if one in- 
creases the Rashba field. Here, the splitting goes from 1 
to 4.2meV if the interface field is increased fourfold. 

Figure 3 illustrates the corresponding level structure 
for two electrons in the QD (full Hamiltonian H + H^e ; 
dashed lines show Hq + H^e, the non-SO case). The re- 
pulsive interaction shifts the ground state upwards by 
~ 5meV, and the exchange shifts the triplet down by 
2meV. Most interestingly, the SO interaction introduces 
a strong mixing of the singlet and triplet transition at 
B ~ 2.7T. The fact that the mixing occurs at relatively 
low field makes that a possibly useful transition for the 
implementation of quantum computing devices. More- 
over, the splitting will also be apparent in the FIR re- 
sponse of QDs, allowing the determination of the various 
SO coupling strengths. 



We have shown that inclusion of all SO terms is es- 
sential in order to obtain a complete picture of the level 
structure in narrow-gap QDs. The combination of strong 
SO couplings and large (and negative) g factor intro- 
duces strong intrinsic mixing of the low excitations for the 
single-particle spectrum. Consequently, the two-particle 
spectrum exhibits strong singlet-triplet coupling at mod- 
erate fields, with significant experimental consequences. 
Observation of FIR mode magnetic dispersion would al- 
low the direct determination of coupling constants. 
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FIG. 3: A: Two particle spectrum vs. B field for full Hamil- 
tonian H + Hee (basis included 190 states, only lowest levels 
shown). With no SO (dotted lines) and at B = 0, ground state 
is a singlet {{L, S} — {0, 0}) at 35meV while first (second) ex- 
cited state is a triplet ({±1, ±1} and {±1, 0}) at 48meV (sin- 
glet ({±1,0}) at 50meV). SO acts against electron-electron 
interaction, as levels are shifted back to energies close to 
non-interacting case. Lowest anticrossing at ~ 2.7T is be- 
tween singlet ground state {0, 0} and lowest excited triplet 
state { — 1, 1}; SO introduces coupling between the singlet and 
triplet states with direct consequences for QD ground state. 
B: Spin Sz for the nine lowest states of the two particle QD. 
Strong mixing induced by SO interaction appears in all states 
~ ST. 
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